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ABSTRACT  (UNCLASSIFIED) 


The  range  prediction  models  REPAS  and  REACT  for  passive  and  active  sonars  respectively,  have  been 
developed  at  FEL-TNO  between  1986  and  1990  .  These  models  and  the  propagation  loss  model  ALMOST 
on  which  REPAS  and  REACT  are  based,  are  documented  in  tins  report. 

The  report  describes  also  the  way  in  which  the  reverberation  is  calculated,  which  is  an  important  part  of  die 
REACT-model.  A  description  is  given  of  two  modules,  v«diich  calculate  reverberation  and  target  echo  level 
as  a  function  of  time  respectively.  Afterwards  time  is  ctmverted  to  detection  range  for  operational  use. 

Finally  a  method  has  also  been  included,  for  the  determination  of  the  threshold  level  for  various  detectitm 
sihiations,  used  in  REPAS  as  well  as  REACT.  With  this  method  also  the  detection  {sobability  can  be 
calculated  versus  detection  range.  Some  examples  of  detection  range  prediction  calculated  by  REPAS  and 
REACT  ate  presented. 
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SAMENVATTING  (ONGERUBRICEERD) 

De  afstandvoorspellingsmodellen  REPAS  en  REACT  voor  achtereenvolgens  passieve  en  actieve  sonars  zijn 
op  FEL-TNO  ontwikkeld,  tussen  1986  en  1990.  Deze  modellen  en  het  propagatieverliesinodel  ALMOST, 
waarop  REPAS  en  REACT  gebaseerd  zijn,  wotden  in  dit  n^tport  gedocumenteeid. 

Het  rapport  beschrijft  ook  de  wijze  waarop  de  reverberade  berekend  is,  die  een  belangtijk  deel  van  bet 
REACT  model  is.  Er  wordt  een  beschiijving  gegeven  van  twee  modules,  die  achtereenvolgens  reverberatie 
en  doels-ecbo-steikte  berekenen  als  functie  van  de  tijd.  Nadeifaand  woidt  voor  operadoneel  gehruik  djd  in 
detecdeafstand  omgezet. 

Tenslotte  is  er  een  mettiode  bijgevoegd,  om  het  drenqtelniveau  te  bepalen  voor  verschillende 
detectiesitoaties,  te  gebniiken  in  zowel  REPAS  als  REACT.  Met  deze  mediode  kan  ook  de  detectiekans 
berekend  worden  versus  detectieafstand.  Er  worden  enkele  vooibeelden  van  detectie-afstandvoot^Uing 
berekend  door  REPAS  en  REACT,  gqrresenteerd. 
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1  INTRODUCnON 

The  aim  of  this  report  is  to  describe  the  following  two  range  prediction  ooodels,  developed 
between  1986  and  1990  at  FEL-TNO: 

REPAS:  Range  Estitnator  for  Passive  sonars  and 

REACT:  Range  Estimator  for  Active  sonars. 

REPAS  and  REACT  both  use  the  relevant  environmental  and  sonar  parameters  in  their 
calculations.  The  ALMOST  propagation  loss  modeL  which  accounts  for  the  environment,  takes  a 
central  place  in  both  nKxlels.  ALMOST  (Acoustic  Loss  Model  for  Operational  Studies  and  Tasks) 
is  the  ptt^agation  loss  model  of  FEL-TNO,  for  operatiotud  use  and  study  purposes,  which  was 
developed  between  1980  and  1988.  The  theoretical  basis  of  ALMOST,  a  modified  ray  dreary,  has 
been  reported  by  Schippers  [1]. 

REPAS  principally  consists  of  an  ALMOST  propagation  loss  calculation,  followed  by  a 
determination  of  the  so-called  Figure  Of  Merite  (FOM),  tthich  is  the  maximum  tolerable 
propagation  loss  for  a  given  passive  detection  situation. 

ALMOST  can  tun  either  with  a  so-called  incoherent  or  a  ciAerent  sound  padi  summation.  The 
incoherent  mode  is  normally  used.  The  coherent  calculation  which  takes  the  phase  of  pressure 
into  account,  generally  gives  an  indication  of  variations  in  propagation  loss  in  case  of  small  band 
detection. 

The  Figure  Of  Merit  (FOM)  is  calculated  m  REPAS  by  means  of  modeb  for  detection  threshold 
and  ambient  noise.  The  ambient  noise  model  uses  an  incoqxnted  data  base.  In  dus  model 
modifications  are  made  on  sea-  and  shipping  noise  widi  die  ALMOST  modeL  to  account  for  sonar 
parameters  and  environmental  conditions.  Also  dte  directivity  index  u  taken  into  account 
REPAS  also  calculates  the  detection  probability  versus  detection  range.  The  detection  module  can 
handle  broad  band  detection  ("energy  detector"  case),  and  fiequency  line  detection  (narrow  band 
frequency  analysu). 

REACT  also  uses  ALMOST,  however  running  die  latter  in  die  "active  sonar"  mode  whidi  has  no 
coherent  alternative.  The  results  of  dus  mode  are  propagation  lots,  arrival  time  and  grazing  angle 
for  each  impmtant  sound  path,  and  for  all  required  ranges.  REACT  contains,  in  coitqnrison  with 
REPAS,  two  extra  modules,  respectively  dealing  widi: 
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the  calculation  of  reverberation  and, 
the  calculation  of  target  echo  level. 

A  data  base  of  sea  surface-  and  bottom  backscattering  is  incoiptvated  in  the  REACT  model.  For 
the  volume  scattering  strength  realistic  default  input  values  are  available  in  the  model,  describing 
the  deep  scattering  layer. 

The  reveiberation  and  echo  contributions  are  calculated  versus  arrival  time,  by  dte  active  scoar 
mode  of  ALMOST.  Arrival  time  is  afterwards  converted  to  the  detection  range  of  active  sonars, 
for  operational  purposes. 

Ambient  noise  is  calculated  in  the  same  way  as  in  REPAS. 

A  detection  module  for  determination  of  detection  threshold  and  detection  probability  for  a  given 
detection  situation,  has  also  been  developed  for  REACT,  which  treat  bodi  CW  and  FM  emissicms. 
Although  REACT  was  principally  developed  for  hull  mounted  sonars,  the  model  could  easily  be 
adt^ted  for  active  adjunct  towed  arrays  and  active  sonar  buoys.  These  various  active  sonar  types, 
each  with  their  own  source  and  receiver  directivity  patterns,  can  be  chosen  by  means  of  a  menu. 

In  the  following  three  chapters  the  structure  and  working  of  ALMOST,  REPAS  and  REACT  will 
be  described. 
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2  ALMOST,  THE  PROPAGATION  LOSS  MODEL 

ALMOST  is  the  propagadon  loss  model  of  FEL-TNO,  which  was  developed,  for  operational  use 
and  study  purposes,  between  1980  and  1988.  The  model  is  based  on  modified  ray  calculations, 
see  ref.  [1]. 

The  model  starts  with  "stylizing"  the  given  Sound  Speed  ^file  (SSP),  which  is  a  frequency 
dependent  data  reduction  process,  by  means  of  which  also  effects  of  diffraction  ate  described.  The 
result  is  a  stylized  SSP  of  which  the  number  of  points  has  been  reduced  to  a  maximum  of  13 
characteristic  ones. 

Then  the  calculation  of  the  propagation  loss  follows.  REPAS  uses  the  propagation  loss  data  per 
range,  by  the  incoherent  or  coherent  calculation  mode.  REACT  uses  propagation  loss,  travel  time 
and  grazing  angle  data  per  ittqtortant  sound  path,  and  per  range  of  the  target  cq.  reverberating  part 
of  the  environment.  For  this  calculation  the  ALMOST  noodel  is  used  in  the  active  sonar  mode. 

In  the  ALMOST  model,  one  can  assume  on  one  hand,  OMNI  reception  in  the  propagation  loss 
calculation,  while  on  the  other  hand  one  can  execute  the  calculation  for  a  vertical  receiver 
sensitivity  pattern,  or  a  vertical  grazing  angle  intaval  of  recepticm. 

Because  the  stylizing  only  requires  input  of  frequency,  SSP  and  bottom  depth,  in  several  cases 
ALMOST  can  be  run  with  the  same  stylizing  ouqnit,  which  saves  computing  time. 


TNOrapoit 


Page 

8 


3  REPAS,  THE  RANGE  PREDICTION  MODEL  FOR  PASSIVE  SONARS 

For  passive  sonars  die  model  REPAS  can  calculate  detection  range  inteival(s)  widi  more  than  a 
given  detection  probability  P^,  or  the  detection  probability  curve  over  a  range  interval  of  interest. 
Id  the  first  case  a  POM  (Figure  Of  Merit)  for  the  given  detection  probability  (Pde^PdeiroM)  ^ 
calculated.  The  detection  module  used  here,  is  the  same  as  to  the  calculation  of  die  detection 
probability  curve. 

We  Stan  with  the  passive  sonar  equation  which  describes  passive  detection  in  general  (see  Urick 
[2]): 


SL-PL-NL  =  DT 


(3.1)  (passive  detection) 


with: 

SL  =  Source  Level  (^level  at  Im  in  receiver  band;  dB  re  |iPa.m) 

PL  s  Propagation  Loss  (dB) 

NL  =  Noise  Level  (self  noise  and  ambient  noise;  in  receiver  band;  dB  re  pPa) 
DT  =  Detection  Threshold  (dB) 


We  observe  in  the  above  equation  that  DT  is  equal  to  the  signal  to  noise  ratio  required  to 
detection. 

In  paragraph  3.1  and  3.2,  the  dependence  of  on  DT  will  be  treated.  Subsequendy  the  passive 
sonar  equation  3.1  will  be  used  to  the  determination  of  the  detection  probability  (Pjei)  versus 
range  curve. 

For  the  prediction  of  detection  range  intervals  widi  a  detection  probability  of  more  than  PeeaoM> 
FOM  can  be  used  (Figure  Of  Merit)  defined  as: 

POM  =  SL  -  NL  -  DTPOM  (3.2) 


with: 

DTPOM  s  Detection  Threshold  DT  to  Pse('*P4c«>M>  vsually  >0.3  (50%) 
We  easily  derive  fiom  eq.  3.1  and  3.2: 
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when: 

PL<  FOM  (3.3) 

In  a  gn^hical  representation  of  eq.  (3.3)  the  so-called  FOM-line  is  drawn  in  the  I^t^>agation  Loss 
(PL)  versus  range  curve.  This  means  that  at  ranges  where  PL  is  smaller  than  die  FOM-value,  a 
detection  probability  of  more  than  PddFOM  expected  (see  fig  1,  where  PdopoM  chosen  =  0.5). 

The  model  REPAS  starts  to  calculate  PL,  the  propagation  loss  versus  range  (curve)  using  the 
ALMOST-model.  Incoherent  or  coherent  sound  path  summation  can  be  chosen,  die  first  being 
used  normally. 

Apart  from  the  omni  directional  calculation,  iiqiut  of  a  vertical  receiver  sensitivity  pattern,  or  a 
vertical  grazing  angle  interval  of  rece{ttion  are  options  in  the  Propagation  Loss  (PL)  calculation. 
FOM  is  calculated  using  the  incorporated  detection  module.  In  general,  this  module  determines 
DT  as  a  function  of  P^„  Pf,  (false  alarm  probability),  the  integration  (processing)  time  T,  and  the 
bandwidth  w  of  the  signal.  This  module  is  used  for  the  calculation  of  POM  by  setting 
Pdei=Pdeif>OM-  Using  the  detection  module  in  another  way  we  can  calculate  P^^  direcdy  finm  DT. 
With  eq.  3. 1  we  then  present  P^  as  a  curve  versus  range,  because  PL  varies  with  range. 

In  the  detection  module  we  distinguish  broad  band  detection,  and  narrow  band  detection,  further 
called  line  detection.  Both  cases  will  be  considered  successively  in  the  next  paragnqrhs.  These 
cases  will  result  in  different  madrematical  expressions  for  P*,.  and  thus  also  for  POM. 

Another  important  term  for  FOM  and  P^  is  NL,  the  noise  level.  This  will  be  treated  in  para  3.3. 
NL  requires  the  horizontal  and  vertical  directivity  pattern  of  the  receiver  (if  present).  This  REPAS 
input  has  been  programmed  menu  driven,  see  para  3.4. 
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3. 1  Broad  band  detection 


Broad  band  signal  detection  is  made  by 
unknown  signal  in  gaussian  noise,  which 
is  determined  by  the  following  equation 
(wT»l;  see  [3],  p.82-88): 


P<kt=  ‘/2erfc( 


X.,-d 


an  energy  detector,  see 
generally  iqiplies  to  titis 
in  an  appronmation  for 


Van  Scbooneveld  [3],  For  the 
passive  detection  situation, 
large  bandwidth-time  products 

(3.4) 


with; 

d  =  10(0  I  DT)  •  V(wT)  (detection  index) 

X  =  ■'/2  crfc-'(2.P(i,)  (sonar  receiver  threshold) 
DT=  SL  -  ML  -  PL 

detection  probability 
Pf,  =  false  alarm  probability 
w  =  bandwidth  of  receiver  (Hz) 

T  » integration  time  of  receiver  (s) 


This  equation  for  P^e,  yields  the  P,|^-curve  versus  detection  range  in  REPAS,  because  PL  is  a 
function  of  range. 

It  is  remarked  here  that  P|^  is  input.  Fck  Gist  detection  it  can  be  chosen  0.0(X)1,  however  after  a 
few  detections  made  of  the  same  target,  it  can  be  increased  up  to  0.(X)1  or  0.01,  to  describe  the 
effect  of  inzooming  by  the  sonar  operator.  This  will  result  in  a  few  dB  lower  DT  and  higher 
For  the  FOM  calculation  of  REPAS,  we  calculate  DTFOM  as  DT  for  die  iiqiut  P^tom 
input  Ph,  by  inversion  of  eq.  3.4.  Then  POM  is  calculated  with  eq.  3.2: 

POM  =  SL  -  NL  -  DTFOM  (3.5) 


with: 


DTFOM=5H»log 


d  -  >/2  [erfc-i(2  Pu>eifc-‘(2  Pi«roM)l 


SL  is  the  source  level,  calculated  by  integration  of  the  spectral  level  over  die  band  of  reception  w. 
In  diis  rqioit  the  noise  NL  is  also  taken  in  die  band  w.  In  fact  the  formulas  3.4  and  3.5  tpply  to  a 
sonar  receWn  with  noise  {newhitenii^.  So  for  SL  we  have  to  integrate  die  level  after 
prewhitening. 
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3.2  Line  detection 

Second,  we  consider  passive  detection  in  a  narrow  band,  of  so-called  line  freqaencies  emitted  by 
targets.  By  definition  the  Source  Level  (SL)  is  here  determined  by  the  anqilimde  of  the  line 
signal.  Just  like  for  broad  band,  this  is  the  level  in  the  band  w. 

The  noise  level  (NL)  is  again  the  integrated  level  over  the  band  w. 

If  a  (reasonably  stable)  line  is  detected  in  a  small  fiequency  band,  we  have  the  case  of  a  (more  or 
less)  known  signal  against  gaussian  noise.  In  general  it  is  expected  that  the  signal  will  fade 
according  to  the  Rayleigh  distribution  for  the  signal  amplitude.  It  is  remariced  here  that  small  band 
filtering  is  in  fact  correlation  with  sines  and  cosines  (Fourier  transform),  so  correlation  with 
tonals.  This  correlation  agrees  with  the  assuiiqition  of  a  completely  known  signal  apart  fi'om 
phase  and  slow  Rayleigh  amplitude  fading  (slow  means  no  fading  during  integration  time).  The 
following  detection  probability  (P^^  holds  for  this  case  (see  van  Schooneveld  [3],  p.74-81): 


^ia~  Pf« 


1/(10“ ‘  •"•wT+D 


(3.6) 


with: 

DT  =  SL  -  NL  -  PL 

This  yields  the  detection  probability  curve  of  REPAS  versus  range,  because  PL  is  a  function  of 
range. 

POM  is  calculated  from  3.2  and  inversion  of  3.6  for  the  input  P^et^PoeiPOM' 


POM  =  SL  -  NL  -  DTPOM  (3.7) 

with: 

DTFOM=  10  'OIog(logPft/IogPfcaoM  -  D '  lOIogwT 
TypicaUy  we  find  POM=10.9dB  for  Pj^-O.OOOl  and  PaaPONrO-5. 
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3.3  Ambient  noise  calculation 

The  noise  level  NL  is  calculated  from  the  spectral  noise.  This  qiectral  noise  consists  of  ambient 
and  self  noise,  both  in  the  receiver  beam  (if  bearofotming  is  ^iplied). 

The  required  spectral  ambient  noise  level  is  calculated  from  an  incoipotated  data  base  with  data 
from  Urick  [2].  The  ambient  noise  data  base  has  the  parameters  wind  speed,  shying  intensity 
and  frequency.  Sea-  and  shipping  noise  are  calculated  separately  for  particular  values  of  these 
parameters  by  interpolation  in  the  noise  data. 

For  a  vertical  receiver  pattern,  the  ambient  noise  is  calculated  as  being  received  in  this  pattern. 
The  ambient  noise  calculation  in  this  case  is  based  on  the  supposition  that  sea  noise  is  generated 
all  over  the  sea  surface,  and  that  shipping  noise  is  generated  by  distant  shipping  beyond  a  distance 
of  27  km.  By  calculating  with  ALMOST,  the  propagation  losses  for  diese  wind-  and  shipping 
generated  sounc*  sources,  for  an  omni-  as  well  as  for  a  vertically  directional  receiver,  both  die  sea- 
and  shipping  noise  level  from  the  data  base  can  be  adtqited  for  the  directional  receiver.  The  result 
is  ambient  noise  due  to  sea  state  and  shipping,  both  adapted  to  environmental  conditions  and 
sonar  parameters,  while  also  the  Directivity  Index  DI  is  taken  into  account. 

Usually  DI  is  calculated.  In  case  of  a  huU  mounted  sonar,  consisting  of  ciicularly  grouped  staves, 
the  receiver  directivity  is  accounted  for  by  subtracting  the  horizontal  Directivity  Index  from  the 
total  calculated  ambient  noise  value.  This  DIj„  is  calculated  by  integration  over  the  horizontal 
receiver  pattern.  For  towed  arrays  however  where  REPAS,  as  an  approximation,  does  not  take  a 
veracal  directivity  pattern  into  account,  the  3  dimensional  DI  is  applied  (^lOlogn,  with  n  = 
number  of  half  wave  lengths  on  receiver  aperture). 

Finally  input  is  required  for  the  sonar  self  noise  in  the  beam.  The  total  spectral  beam  noise  level  is 
subsequently  calculated  as  the  intensity  sum  of  ad^ned  sea  noise,  adapted  shying  noise,  both 
minus  DI,  and  self  noise  in  the  beam. 
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3.4  Menu  driven  choice  of  sonar  type 

The  passive  sonar  input,  consisting  of  the  vertical  angle  interval  or  receiver  directivity  pattern, 
and  the  horizontal  pattern,  is  menu  driven  for  easier  use.  In  particular  the  different  calculation 
methods  of  these  patterns  and  their  application  in  REPAS  for  either  towed  array  or  HMS,  require 
such  menu  driven  input.  The  sonar  menu  has  the  following  form,  where  the  menu  number  can  be 
chosen  in  REPAS; 

code: 


1, 

Towed  receiver  array 

pta 

2. 

Single  hydrophone 

shy 

3. 

Vertical  receiver  array  (for  instance:  passive  buoy) 

via 

4. 

Hull  mounted  sonar 

hms 

See  for  more  details  about  the  calculation  of  patterns  for  linear  arrays  and  HMS:  paragraph  4.4. 

As  indicated  in  the  menu  list  REPAS  can  be  used  for  towed  array  sonar,  buoys,  VDS  and  HMS. 

A  REPAS  plot  is  shown  in  fig  1  for  the  passive  mode  of  a  hull  mounted  sonar,  calculated  for  the 
sound  speed  profile  of  fig  2,  a  deep  winter  profile.  The  heading  shows  the  various  REPAS  (and 
ALMOST)  parameters,  starting  with  the  menu  choice.  The  plot  shows  the  propagation  loss  curve 
(PL)  with  FOM  line  on  one  hand,  and  on  the  other  band  at  the  bottom  the  curve,  versus  the 
passive  detection  range.  It  is  remarked  that  accarding  to  the  definition  of  FOM,  P^  equals 
^deiFOM  to  fig  1)  at  the  range  where  PL  equals  FOM. 

In  addition  the  main  contributing  ALMOST  propagation  modes  (see  ref  [1])  are  plotted  at  8 
different  ranges  on  the  range  interval.  Further  the  total  number  of  modes  for  this  situation  is 
presented  at  the  top  of  the  cadre.  We  see  in  fig  1  that,  for  the  winter  profile  of  fig  2,  surface  duct 
propagation  is  dominant. 
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4  REACT,  THE  RANGE  PREDICTION  MODEL  FOR  ACTIVE  SONARS 

The  model  REACT  calculates  for  active  sonars,  detection  range  interval(s)  with  more  than  50% 
detection  probability  (Pdet),  or  the  P^-curve  over  a  range  interval  of  interest.  For  the 
determination  of  the  50%  detection  range  the  Minimum  Detectable  Signal  (MDS)-cutve  is 
calculated  for  P,„=0.5  (50%).  This  ciuve  is  compared  with  the  calculated  target  echo  level  curve. 
If  the  echo  level  is  larger  than  MDS  the  detection  probability  is  more  than  50%. 

Active  detection  is  described  by  the  active  sonar  equation  (see  Urick  [2]): 

SL  -  PL,„  +  TS  -  PLj^  -  BGL  =  DT  (4.1) 

^  target  uiget  ^ 

V 

TEL 


with: 

SL  =  Source  Level  of  sonar 

TS  =  Target  Strength  (=''nominal  secondary"  SL  of  target) 

BGL  =  Background  level  (noise  and  reverboaticm) 

DT  =  Detection  Threshold 
TEL  =  Target  Echo  Level 

We  observe  in  the  above  equation  that  DT  equals  the  signal  to  noise  ratio  required  for  detection. 

In  para  4. 1 ,  the  dependence  of  P^jat  on  DT  will  be  treated.  Subsequently  the  active  sonar  equation 
4. 1  will  be  used  for  the  determination  of  the  detection  probability  (Pget)  versus  range  curve. 

For  the  prediction  of  detection  range  intervals  with  mote  than  50%  detection  probability  MDS  can 
be  used  defined  as: 

MDS  =  BGL  +  DT50  (4.2) 


Like  DTFOM  in  REPAS,  ch.  3,  DT50  is  DT  for  Paa=0.5. 
We  easily  derive  fi'om  the  eq.  4. 1  and  4.2: 

Pja,>0.5(50%) 

when; 


TEL  >  MDS 


(4.3) 
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Both  the  TEL-  and  MDS-ciuve  can  be  plotted  by  REACTS  plot  routine,  from  which  the  interval 
of  more  than  50%  detection  probability  can  be  read  out  (see  fig  3).  DTSO  will  be  calculated  in 
para  4.1. 

In  eq.  4.1  we  see  propagation  loss  (PL)  to  and  from  the  target.  In  the  same  way  there  is  PL  to  and 
fi'om  any  reverberating  part  of  the  environment,  cmsisting  of  sea  surface,  volume  (Deep 
Scattering  Layer  DSL)  and  bottom.  For  DSL  the  input  is  9  layers.  PL  is  calculated  for  these 
mechanisms,  and  for  the  two  ways  separately,  ^tplying  ALMOST.  REACT  qrplies  ALMOST  in 
the  so-called  active  sonar  mode.  The  stylizing  of  ALMOST  is  used  in  the  same  way  as  for 
REPAS,  but  in  REACT  always  for  one  frequency  per  calculation.  The  ALMOST  propagation  loss 
calculation  yields  here  per  important  propagation  padi:  PL,  arrival  tiine,  and  grazing  angle  at 
target,  for  all  target  ranges  cq.  ranges  for  the  reverberation  (RL)  calculation.  For  respectively 
target  and  environment  (reverberation),  all  combinations  of  paths  "to”  and  "fiom"  are  made,  in 
respectively  the  TEL-  and  the  RL  calculation.  REACT  calculates  "quasi  monostatic",  which 
means:  source  and  receiver  at  same  geographic  position  with  option  of  different  depths.  Also  the 
vertical  patterns  of  reception  and  emission  may  be  different.  In  ctmtrast  to  REPAS,  a  simple 
vertical  angle  interval  for  source  or  receiver  is  not  allowed  in  REACT  because  the  prediction  may 
become  poor  if  the  effects  of  sidelobes  are  neglected,  which  will  be  shown  further  in  this  chapter, 
in  para  4.5.  The  calculations  are  made  versus  the  arrival  time  which  is  the  sum  over  a  path  "to” 
and  a  path  "from".  In  general  there  is  backscatteiing  with  different  angles  of  incidence  and 
scattering.  In  particular  for  the  reverberation  calculation  this  has  to  be  taken  into  account  Further, 
there  is  for  REACT  no  alternative  coherent  calculation  mode,  since  cdierence  will  be  lost  by  die 
scattering  mechanism. 

A  few  remarks  still  have  to  be  made  about  the  reverberation  calculation.  For  surface  scattering, 
the  distributed  secondary  sound  sources  are  placed  at  a  dq>th  of  0.1m,  because  ALMOST  does 
not  accept  an  input  of  Om.  (In  general  die  acoustic  pressure  will  be  zero  at  the  sea  suifiux.)  The 
so-called  surface  decoupling  loss  has  been  omitted  in  die  active  sonar  mode  of  ALMOST  to  get  a 
realistic  description.  Moreover  as  there  are  allways  two  nearly  congruent  rays;  with  and  widiout 
reflection,  the  (linear)  scattering  strength  has  been  divided  by  2,  which  is  also  of  tpplicadoo  for 
the  bottom  scattering  data. 

Finally  there  are  some  renuuks  about  die  effective  surface-  cq.  volume  part  to  udiich  the 
scattering  data  are  ipplied.  Although  reverberation  is,  in  fact,  die  convolution  of  distribalBd 
scatterers  with  the  emitted  pulse  shape,  this  conplicated  calculadtm  is  qiinoximated  by  taking 
parts  of  size  cT/2  in  the  direction  of  the  np  paths.  However,  die  angle  of  bieidence  and  die 
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division  of  the  environment  in  range  intervals,  also  detennine  the  horizontal  size  of  die  scattering 
part,  which  is  accounted  for  in  the  reveiberation  calciilaticHi. 

The  single  paths  fcom  source  to  receiver,  SRP-  and  BSP  mode  propagation  in  ALMOST  terms 
(see  ref.  [1]),  are  calculated  versus  time,  by  application  of  the  active  smiar  mode  of  ALMOST. 
The  well-known  "fathometer  effect"  is  included  in  this  type  of  detectimi  background. 

These  single  paths  are,  in  many  sonars,  received  in  sidelobes,  which  results  in  suppression  of  the 
effect.  For  towed  receiver  arrays  the  REACT  approximation  does  not  apply  a  vertical  receiver 
pattern.  However  for  steering  angles  which  are  not  broadside  die  sidelobe  level  in  die  vertical 
direction  is  accounted  for. 

For  a  single  target,  one  emitted  pulse  can  result  in  several  echoes  splitted  in  time.  REACT 
supposes  a  number  of  identical  targets  over  an  input  range  interval  (maximum  1(X)  targets).  TEL 
is  calculated  by  taking  at  each  time  bin,  the  strongest  of  the  echo  levels  of  all  these  targets. 

The  reverberation  is'  calculated  for  the  same  range  interval.  This  part  of  the  environment  (sea 
surface,  -volume  and  -bottom)  is  divided  in  a  number  of  range  intervals  (max.  100)  each 
producing  (back)scattered  intensity.  This  number  of  ranges  is  taken  equal  to  the  number  of  target 
ranges. 

Afterwards  the  arrival  time  scale  is  replaced  by  a  conesponding  detection  range  scale  in  the 
operational  plot. 

The  remainder  of  this  paragraph  will  concentrate  on  modeling  of  detection  threshold  and  doppler, 
followed  by  general  remarks  about  the  backscatteiing  data  base  used  for  die  reverberation 
modeling.  After  this  whole  basic  description  of  REACT,  the  menu  driven  iiqnit  for  diffeient  sonar 
types  will  be  described,  and  some  exatiqiles  of  REACT  predictions  will  be  given. 


TNO  report 


17 


4.1  Detection  threshold 

The  modeling  of  Detection  Threshold  DT  depends  on  die  statistics  of  echo  signal  and 
background,  and  the  type  of  sonar  receiver.  This  problem  can  also  be  fcimulated  as  calculating 
the  detection  probability  as  a  function  of  DT,  which  is  in  fact  the  signal  to  noise  ratio. 

In  a  modem  active  sonar  the  echo  signal  is  correlated  with  die  known  emitted  signal  shtqie.  This 
is  the  detection  case  of  a  completely  known  signal  against  a  background  of  noise  and 
levedieration.  The  echo  signal  will  possess  fading.  This  is  mathematically  the  same  detection  case 
as  for  passive  line  detection  (para  3.2). 

To  calculate  DT  for  the  background  of  noise  plus  Kverberadon,  first  the  two  limit  cases  will  be 
considered:  background  of  noise  only,  and  of  reverberation  only.  Finally  the  general  mixed 
background  will  be  considered. 

The  formula  of  Pj^,  the  detection  probability  for  the  case  of  a  known  signal  with  slow  Rayleigh 
fading  in  gaussian  noise  is  as  follows  (see  Van  Schooneveld  [3],  p.74-Sl,  and  also  eq.  3.6): 


■'del*  •'ft 


l/(10“‘“^+l) 


(4.4) 


with: 

DT=  TEL  -  NL+10*‘«log(wT) 

NL=  noise  level  in  band  w 
w  =  bandwidth  of  receiver  (Hz) 

T  =  integration  time  of  receiver  (s)  =  pulse  length 

The  only  difference  with  eq.  3.6  is  that  the  wT  product  is  larger  than  1  here.  NL  is  calculated,  like 
in  REPAS,  accounting  for  sonar  self  noise  as  well  as  ambient  noise  level  in  the  vertical  receiver 
pattern  (input:  wind  ^leed,  shaping  intensity;  see  para  3.3). 

Because  the  ambient  noise  level  is  calculated  as  dw  omni  directioiial  level  received  in  the  vertical 
receiver  pattern  of  the  sonar  (see  REPAS  para  3.2),  this  value  is  snbseqnentfy  decreased  with  the 
horizontal  Directivity  Index  DI.  The  latter  is  calculated  by  integration  over  the  horizoatal  receiver 
pattern.  This  method  is  followed  for  hull  mounted  sonars,  consistiDg  of  staves,  cncularly  groi^ied, 
and  for  active  sons  buoys.  However  an  exceptioa  is  made  for  die  towed  receiver  array  (with 
active  adjunct).  This  has  no  vertical  paMem  in  the  REACT  ^iproximation,  so  DI  is  calculated  here 
based  on  the  number  of  half  wave  lengths  n  on  the  atn^  length  (DIwKNogn). 
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Subsequently  DI  is  subtracted  from  sea-  and  shipping  noise  level,  and  then  conibined  widi  die 
sonar  beam  self  noise. 

If  the  above  known  signal  must  be  detected  in  a  background  of  reverberation  (RL),  possessing  a 
Rayleigh  distribution,  we  calculate  in  a  good  approximation  the  same  farmula  for  however 
without  background  reduction  of  wT  like  gaussian  noise  in  4.4; 


p  _  '/(10“‘'>T+i) 

P M-  Pfl 


(4.5) 


with: 

DT=TEL-RL 
RL=  reverberation  level 


The  matched  filter  receiver  causes  time  conqnession  of  the  emitted  pulse  signal  from  T  down  to  a 
time  1/w  (see  ref.  [3]).  This  results  in  a  reduction  of  RL  with  lOlogwT. 

For  the  general  detection  case  of  a  background  of  noise  plus  leveiberation  the  following 
formulation  for  P^,  is  adopted,  in  agreement  with  the  above  limiting  cases  (eq.  4.4  and  4  J); 


P*,=  Pf 


l/(lop.ior+i) 


(4.6) 


with; 

DT  =TEL-BGL 

BGL  =  intensity  summation  of  RL  and  NL-lOlogwT 

In  the  following  this  expression  for  will  be  ^iplied  to  the  two  main  pulse  types  for  active, 
matched  filter  sonars: 

FM;  nequeucy  Modulated  pulse  and 
CW:  Continuous  Wave  pulse. 

In  the  latter  case  (CW)  w  is  relatively  smaB,  wwl/T.  So  in  the  CW  case,  die  product  wT  in  eq.  4.6 
equals  1.  Generally,  a  doppler  shift  in  frequency,  of  the  received  signal  is  caused  by  sonar  and 
target  speed.  For  each  frequency  interval  w  of  the  whole  leceptioo  band  for  dopider  shifted 
echoes,  a  sepat^  matdied  filter  is  needed.  This  receiver  type  is  usually  jndicMed  as  a  'doppler 
filter  bank"  receiver.  It  makes  in  frKt  the  Pouirier  nnsformatioD  of  die  signal  in  die  entire 
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receiver  band.  Of  course  there  is  no  titne  concession  in  die  CW  case.  So,  compand  with  FM,  die 
resolution  of  the  active  detection  ranges  found  widi  CW,  is  usually  much  worse. 

It  must  be  remariced  that  also  in  FM  sonar  receivers  often  a  few  doppter  shifted  replicas  are  used 
for  the  optimum  matched  filter  ouqiut-  REACT  supposes  such  optimized  receivers. 

The  MDS  curve  is  calculated  from  DTSO,  see  eq.  4.2.  DTSO  follows  from  inversion  of  eq.  4.6, 
like  3.6,  of  course  for  both  CW  and  FM  (CW  is  in  fact  FM  with  wT*l): 


DTSO  =  10  '0log0ogP,^ogP,fc,  -  1)  (4.7) 

with: 

Pto  =  0.5 

Typically  we  find  DT5(^10.9dB  for  Pfa^.OOOl. 

There  are  however  older  CW  sonars  which  do  not  use  a  doppler  filter  bank.  The  recqidon  band 
for  the  CW  mode  is  designed  such  wide  that  usuaUy  all  possible  doppler  shifted  target  echoes  are 
received.  This  relatively  simple  detector  type  can  be  described  as  the  energy  detector  (see  [3], 
p.82-88,  and  also  eq.  3.4): 

Part-‘/2«fc(^^j 

with: 

d  a  lOtO  l  DT) 

X  »V2erfc-K2Pa) 

DT  -TEL-BGL 
BGL  >  intensity  sum  of  RL  (wifiiout  time  ctmopression)  and  NL  -  S*%gwT 

In  REACT  this  type  of  sonar  has  been  called  "CW-enogy”  detectian.  MDS  follows  in  this  case 
again  from  4.2,  with  DTSO  calculated  by  inversion  of  4.8  (see  also  3.5): 

DT50«  lO  'OlogX.  (4.9) 

A.-V2-erfc-i(2-Pfc) 


(4.8) 

(detection  index) 

(sonar  receiver  threshold) 


where  Paei^.S  has  already  been  substituted. 
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4.2  Doppler  effects 

It  is  possible,  that  a  target  echo  possesses  a  doppler-shifi  (m  the  ffequeacy  scale  cooqMied  to  the 
reveiberatitHi,  due  to  the  target  qteed.  ff  target  echo  and  reverheratioa  do  not  fall  in  die  same 
doppler-filter  band  (width  w),  detection  has  to  be  made  against  a  noise  background  NL  only.  In 
particular  for  CW,  wifh  small  w  values,  diis  results  in  a  more  favourable  detection  situation. 

For  shifts  smaller  than  w,  a  decreased  reverberatimi  level  is  taken  into  account  based  cm 
interpolation. 

Also  the  speed  of  the  sonar  vessel  is  of  in^xirtance  for  RL,  causing  the  following  conqilicated 
effect.  Reverberation  contributons  from  different  horizontal  directions  possess  different  doppler 
shifts.  Much  of  this  reverberadon  is  received  in  sidelobes.  Some  reverberation,  mostly  frmn 
sidelobes,  may  fall  in  the  doppler  band  of  the  main  target  echo  Cmain”  in  case  of  echo  splitting). 
This  effect  has  been  modeled  in  REACT,  by  the  calculation  of  a  horizontal  angle  interval  for 
which  the  reverberation  has  a  doppler  shift  conqiarable  to  that  of  the  target  echo.  The  average 
direction  together  with  the  effective  width  of  this  horizontal  angle  interval  are  calculated.  The 
reverberation  contribution  is  proportioiial  to  this  width. 

In  REACT,  doppler  shifts  of  reverberation  and  single  padis  are  calculated  relative  to  the 
fi-equency  shift  of  the  main  target  echo.  The  latter  will  in  general  be  shifted  in  the  absolute  sense 
(Af)  if  there  is  a  target  speed  component  in  the  main  beam  (Urget)  directica.  REACTT  calculates 
the  doppler-shift  relative  to  the  target  echo  Af^  fiom  the  absolute  doppler  shift  Af,  for  each 
reverberation  contribution. 

The  absolute  doppler  shift  of  the  target  echo  is  determined  by  the  projection  of  sonar  and  target 
speed  on  the  vertically  inclined  sound  paths.  Fbr  the  absolute  doppler  shift  of  reverberation  the 
"target  speed"  must  be  taken  zero: 

^  •  «o»  ®  *  V,,)  •  (cos  0„-K»s  9^ /c  (4.10) 

with; 

Af  >  Doppler  fiequency  shift 
f  sftequeocy 

s  speed  of  sonar  ship 

»  speed  of  target  in  target  direction  (bearing;  >0>leaving; -0  ftir  RL) 

O  s  bearing  angle  of  target  (hmizontal) 
c  a>soundspeed 

9,0  -  grazing  angle  of  path  to  target  (vertical) 

Oftm  *  grazing  angle  of  padi  fiom  target 
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The  single  paths  (fathometer  effect)  possess  no  absolute  doppler  shift,  as  source  and  receiver  have 
the  same  speed,  and  the  reflecting  horizontal  bottom  is  right  below  source  and  receiver. 

It  should  finally  be  rematked  that  for  moving  targets  no  or  very  little  advantage  can  be  achieved 
by  using  doppler-shiiis  in  case  of  FM-pulses.  This  is  due  to  the  relatively  large  bandwiddi  w  of 
the  firequency  sweep  of  such  pulses. 
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4.3  Scattering  data 

REACT  requires  scattering  data  for  dte  surfaces  (sea  surface  cq.  -bottom)  and  the  volume.  For  the 
surfaces  the  usually  different  angles  of  incidence  and  scattering  in  REACT,  as  we  have  seen 
above,  are  parameters  for  the  scattering  stiengdi.  FOr  surfaces  the  target  stroigdi  per  m^  is 
required,  and  for  the  volume  the  target  strength  per  m^,  the  latter  independent  of  the  angles. 

Such  detailed  data  for  the  surfaces  are  not  available  in  the  literature.  It  should  be  noted  that 
bottom  type  and  frequency  are  also  parameters  for  bottom  scattering.  In  the  same  way  wind  speed 
and  frequency  are  also  parameters  for  sea  surface  scattering. 

In  the  literature  only  data  for  real  backscattering  are  available,  in  stead  of  scattering  into  various 
angles.  So  these  data  have  to  be  adapted  for  use  in  REACT,  to  handle  the  different  angles  of 
incidence  and  scattering.  This  adaptation  method  will  be  treated  below.  First  however  some 
remarks  about  the  data  from  the  literature. 

A  database  for  sea  surface-  and  bottom  backscattering  strength  as  a  ftinction  of  grazing  angle  and 
frequency  has  been  constructed.  In  fact  die  sea  surface-  and  bottom  backscattering  curves  (versus 
angle)  have  been  read  out  in  their  charac  teristic  points,  and  added  to  the  data  base. 

Sea  surface  data  have  been  taken  from  Urick  [2].  Arrays  of  backscattering  curves  versus  angle  are 
given  with  wind  speed  as  parameter.  The  observed  frequency  dependence  in  these  data  is  not  vety 
large.  Therefore,  a  data  set  representing  these  curves  is  conqtosed  for  each  of  the  following  three 
adjacent  frequency  bands: 

1)  0.02-  2.5  kHz 

2)  2.5  -  13  kHz 

3)  13  -  31  kHz, 

thus  covering  a  large  frequency  range  of  iqiplication:  20  Hz  to  31  kHz.  This  means  that  no 
frequency  intopolaticHi  is  made  here  in  REACT. 

Bottom  backscattering  curves  for  bottoms  widi  diffnem  roughness  (-type),  are  gives  in  Utkk  [2] 
in  a  similar  way.  The  measure  of  bottom  roughness  has  been  linked  to  the  MOS-nutnber,  already 
used  in  ALMOST  (see  [1]),  which  could  only  be  done  in  a  global  qtproximate  way  (die  MGS 
number  is  principally  used  to  describe  die  forward  bottom  reflectkm  loss).  The  obtained  ain^  of 
backscattering  curves  versus  angle  widi  die  MOS-number  as  die  adopted  parameter,  are  used  in  a 
similar  way  as  the  sea  surface  backscattering  data.  The  frequency  intervals  of  qiplkation  of  these 
bottom  data  are  also  large,  so  for  practical  reasons  the  same  frequency  bands  have  been  diosen  as 
for  die  sea  surface. 


TNO  report 


Page 

23 


The  calculation  of  the  backscattering  for  the  given  wind  speed  cq.  MGS-number,  from  the  above 
constructed  data  base  is  perfonned  by  calculation  of  an  inteipolated  curve  (versus  angle)  for  the 
required  parameter  values.  The  data  sub  set  for  the  relevant  frequency  band  is  qrplied  here. 

The  backscattering  data  (for  equal  angles  of  incidence  and  scattering)  are  adiqtted  in  REACT  to 
detennine  the  scattering  in  any  direction,  which  is  required  for  the  calculation  of  die  reverberation 
level.  This  is  achieved  by  making  use  of  Lambert's  scattering  law  for  angular  intensity 
distribution  (see  Utick  [2]); 

sin  6. 

ss(0i,e.)  =  ss(ei)  — ^  (4.ii) 

sin  6j 


with: 

SS  =  Scattering  Strength  of  sea  surface  cq  -bottom,  per  m^ 

Oi  =  incident  grazing  angle  (vertical) 

6,  =  angle  of  scattering 

The  database  can  easily  be  updated  if  better  data  become  available.  Because  of  lack  of  scattering 
strength  data  for  small  angles,  extrapolations  have  been  made  according  to  the  behaviour  of  high 
frequency  data  in  ref  [2].  With  these  data  REACTS  results  conq»re  well  with  measurements  in 
literature. 

The  volume  backscattering  data,  given  as  a  function  of  depth,  form  a  default  irqrut  for  a  Deqi 
Scattering  Layer  (DSL).  The  volume  scattering  values  ate  only  9  defrnilt  figures  for  the  9  sub 
layers  in  REACT,  which  can  be  changed  if  requited. 

The  data  taken  from  Utick  [2],  ate  princqwUy  independent  of  grazing  angle.  A  frequency 
dependent  trend  of  the  mean  scattering  level  of  the  DSL,  as  given  by  Urick  [2],  is  also 
hrqrlemented  in  REACT.  This  trend  forms  a  good  defrmlt  descrqrtion  of  the  DSL,  but  it  can 
sinqrly  be  deleied  for  qiecial  cases. 

The  depth,  thickness  and  scattering  level  of  the  DSL  can  easily  be  aibqMed  in  REACT  by 
changing  the  sub  layers  of  constant  scattering  strengdi;  Fbr  instance  a  depfii  change  of  the  whole 
DSL  may  be  important  for  a  better  descrqrtioo  of  the  day  and  night  situation. 
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4.4  Menu  driven  choice  of  sonar  type 

REACT,  at  this  stage,  can  be  run  fts-  HMS,  active  adjunct  and  active  sonar  buoys.  It  is  of  great 
in^rtance  that  the  input  source  and  receiver  directivity  patterns  are  chosen  in  die  ri^t  way. 
There  are  two  possibilities,  calculation  of  the  patterns  firoin  the  geometry  of  (in  most  cases) 
distributed  projectors  and  hydrophones,  or  reading  the  patterns  fixim  a  file.  The  menu  shown  next, 
particularly  assists  when  the  option  is  chosen  of  calculating  the  patterns.  The  menu  also  plays  a 
role  in  the  calculation  of  the  Directivity  Index  DI  and  the  sidelobe  level  at  which  the  single  paths 
(fathometer-effect)  are  received.  This  means  that  the  input  menu  number  is  used  in  all  cases  of 
REACT  calculations. 

The  calculation  of  the  patterns  is  done  for  linear  source  and  receiver  arrays  in  an  exact  way, 
however  supposing  no  shading. 

HMS  sonars  have  vertical  staves  which  are  grouped  in  a  (horizontal)  circle.  Supposing  that  a  half 
circle  of  staves  emits  cq.  receives,  the  horizontal  pattern  is  calculated  by  using  a  linear  rqtlacing 
array  coinciding  with  the  diameter  of  the  circle,  and  filled  with  elements  at  a  selectable  spacing 
("input  spacing"),  which  is  usually  taken  as  half  a  wave  length.  The  beam  sluqte  is  always 
calculated  for  the  broadside  case.  This  linear  replacing  array  yields  in  broadside  two  symmetric 
beams  of  which  one  is  deleted  in  the  calculation,  because  such  a  second  beam  is  not  present  in  the 
HMS  beam  forming.  In  order  to  calculate  realistic  patterns  a  selectable  minimum  sidelobe  level  is 
taken  into  account  in  all  pattern  calculations  used  in  REACT  (and  REPAS). 

The  vertical  HMS  pattern  is  calculated  starting  with  die  pattern  of  one  stave.  Linearly  grotqied 
elements  ate  assumed,  at  equal  "input  spacing",  in  agreement  with  die  geometry  of  a  HMS 
(stave).  Subsequendy  this  pattern  is  mult^lied  by  the  vertical  patlem  of  die  (emitting  cq. 
receiving)  horizontal  half  circle.  This  second  pattern  is  again  calculated  nsmg  a  replacing  linear 
array.  This  linear  airay  is  chosen  half  the  length  of  the  circle  diameter,  horizontally,  in  the  beam 
direction,  with  the  elements  at  die  selected  "isput  qwcing",  and  steered  in  endfire  direction. 

In  case  diat  die  sonar  does  not  use  exacdy  half  die  number  of  staves  for  die  beamforming,  die 
sipposed  acoustic  diameter  can  be  chosen  somewhat  smalkr.  The  in  dns  way  caknlated  HMS 
patterns  appear  to  ag^  very  well  widi  data  of  the  mamifactnrer.  The  advantage  of  pattern 
calculation  is  that  fiequency  changes  are  automatically  taken  into  account,  given  the  geometry. 
The  sonar  menu  kxdcs  as  follows,  tdiere  the  menu  number  can  be  chosen  in  REACT: 


code; 
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1.  Towed  receiver  array/ vertical  source  array  atd 

(for  instance;  ATAS,  active  adjunct) 

2.  Towed  receiver  array/  single  source  aim 

3.  Vertical  receiver  array/ vertical  source  array  bvd 

(for  instance:  active  buoy) 

4.  Vertical  receiver  array/  single  source  bvm 

5.  HMS/  Omni-  cq.  TRDT  emission  hms 

6.  HMS/  single  beam  emission  hmb 


It  is  remarked  that  for  a  HMS  used  in  the  so-called  TRDT  mode,  the  prediction  has  to  be  made 
using  the  same  menu  (hms)  as  for  OMNI,  however  with  the  appropriate  higher  source  level. 

It  is  possible  to  expand  the  menu  if  required. 


I 


TNO  report 


Page 

26 


4.5  Examples  of  REACT  calculadons 

Some  examples  of  REACT  plots  are  shown  in  fig  3  to  8,  for  an  active  adjunct  low  fiequency 
sonar,  an  ATAS  sonar  (Active  Towed  Array  Sonar),  and  a  hull  mounted  sonar  for  which  the 
REACT  model  was  originally  developed.  The  Sound  Speed  Profile  (SSP)  given  in  fig  2,  is  a 
typical  deep  winter  SSP.  Sonar  and  target  are  situated  within  the  90m  deep  Surface  Duct.  As 
indicated  in  the  heading  of  the  plots  the  full  array  geometries  of  the  sonars  chosen  have  been 
taken  into  account. 

Successively,  the  sonar  menu  choice  and  all  REACT  parameters  of  sonar  and  environment,  are 
given  in  the  heading  of  the  REACT  plots.  The  curves  are  plotted  versus  range,  and  echo  levels  are 
compared  with  MDS.  All  plots  are  of  the  so-called  operational  type  except  for  fig  5  which  is  of 
the  so-called  theoretical  type.  Here  the  echo  splitting  in  time  is  presented  explicitely  using 
different  symbols  for  the  target  ranges.  The  maximum  number  of  target  ranges  is  only  20,  in  stead 
of  100  for  the  operational  plots.  In  the  theoretical  plot  the  reverberation  level  is  plotted  separately, 
also  indicating  the  main  contributing  reverberation  mechanism(s): 

S  =  surface; 

B  =  bottom 

L  =DSL. 

For  the  operational  plots,  P^  is  presented  as  a  function  of  range  as  well.  It  is  remarked  that 
equals  0.5  at  the  range  where  echo  level  and  MDS  ate  equal. 

The  importance  of  realistic  beam  patterns  in  stead  of  approximations  by  rectangular  beam  shapes 
is  illustrated  by  fig  7  and  8.  In  fig  8  a  rectangular  horizontal  beam  sh^te  without  sidelobes  is 
used,  with  a  realistic  width  of  11  degrees  for  the  HMS  sonar.  Conqrated  to  fig  7,  using  fire 
complete  horizontal  pattem,  a  big  difference  in  reverberation  level  is  shown,  even  for  the  hi^ 
target  speed  of  6m/s.  A  positive  speed  value  means  that  the  target  moves  away  fiom  the  sonar 
ship,  in  the  beam  direction,  with  the  given  speed,  while  a  negative  value  means  tqrproaching.  In 
fig  7,  doppler  shifted  reverberation  is  received  in  sidelobes  only,  but  diis  level  appears  to  be  of 
importance  in  the  REACT  prediction.  In  fig  8  this  level  is  not  present  due  to  the  rectangular  beam 
shqre  with  no  sidelobes,  resulting  in  a  too  optimistic  prediction.  It  is  remarked  that  a  change  in 
Target  Strength  TS  sinqrly  results  in  a  vertical  shift  of  the  Target  Echo  Level  curve. 

In  fig  10  to  13  die  REACT  calculatimis  of  fig  S  to  8  have  been  rqpeated  for  a  summer  case,  for 
which  the  SSP  is  given  in  fig  8.  We  see  that  the  high  reverberation  contribution  fm  die  winter 
case  (fig  S  to  8),  due  to  surface  reverberation,  is  not  present  in  die  summer  case.  Also  die  echo 


* 


TNO  report 


P»ge 

27 


levels  are  much  lower  due  to  bad  propagation  conditions  (no  surface  duct).  The  small  maximum 
in  fig  11  at  about  3km  represents  bottom  echoes  from  very  close  targets  as  shown  in  the 
theoretical  plot  (fig  10).  In  fig  12  and  13  this  feature  is  not  present  anymore  because  at  the  target 
speed  of  6m/s  there  is  a  large  doppler  shift  difference  between  direct  and  bottom  reflected  paths, 
due  to  the  steep  grazing  angles  of  the  latter.  This  is  contrary  to  the  case  with  target  speed  Om/s  of 
fig  10  and  1 1 ,  where  the  doppler  shift  diffeience  between  direct  and  bottom  paths  is  such  that  the 
echoes  fall  partly  in  the  same  doppler  band,  causing  the  small  noaximum  at  about  Skm. 
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5  CONCLUSIONS 

A  complete  description  of  the  passive  and  active  sonar  range  prediction  models  REPAS  and 
REACT  has  been  presented  in  this  leport,  together  with  some  exanqtles,  for  HMS  and  active 
adjunct  sonars,  in  two  different  seasons. 

The  REPAS  model  uses  both  propagation  loss,  determined  by  the  already  tested  ALMOST  model, 
and  a  detection  module,  to  calculate  the  final  passive  stmar  range  prediction. 

The  REACT  model  also  uses  the  detection  module,  but  calculates  in  addition  the  reverberation 
level  with  the  aid  of  a  special  calculation  mode  of  ALMOST.  ALMOST  is  used  here  for  the 
calculation  of  propagation  losses  and  arrival  times  for  separate  sound  paths.  The  method  of 
calculation  is  also  applied  to  the  calculation  of  the  splitted  target  echo  level.  In  addition  doppler 
effects  in  all  sound  paths  are  taken  into  account  for  the  reverberation  as  well  as  for  the  echo.  An 
incorporated  data  base  of  environtnetttal  backscattering  data  is  used  for  the  reverberation 
calculation. 

Both  models  further  make  use  of  an  ambient  noise  model,  using  the  receiver  beam  pattern  and  an 
incorporated  ambient  noise  data  base. 

The  menu  driven  input  for  both  REPAS  and  REACT  allows  for  operational  use  at  sea  arxl 
facilitates  their  use  in  studies.  The  menu  takes  care  of  the  specific  sonar  input,  which  is  rather 
complicated,  due  to  the  directivity  patterns  resulting  from  the  sonar  geometry.  In  this  way  the 
input  of  false  or  utuealistic  parameters  is  avoided. 
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Fig.  2:  Historical  sound  speed  prorile  of  winter  type. 


P4t  UrtI  TARGET  ECHO  LEVEL  .od  HOS  50  pel  ■  i^) 


TNOnpoct 


Page 

32 


INPUT  PARAMETERS 


SONAR/SSP 

•id 

aph^AA^I  202 

FM  pyUa  > 

BOOO  ai«/60  H: 

* 

750. 

.0  H« 

ncs  Ppov. 

■  2 

SRC.P.i.Vt 

ALf 

1-00*.  Xe.OOOa 

Uata^  DpiK  1 

5200  .. 

no  REC • Pat  • 

.  V 

WIa^  Spaad  > 

7-0  */. 

REC.P.t.Hi 

Atr 

).00*.92k.t30« 

Sk1pp>  ini • 1 

1 .0  (  LOU 

SRC • Dpih  * 

70  » 

SL 

■  218  dB 

REC'Opih  t 

o 

f 

TS 

1  10  dB 

TGT • Opih  * 

60  • 

Spd.SON/TGT 

«  6.0  /  -3.0 

Hg.  3: 


REACT,  using  ALMOST;  towed  amy  with  active  adjunct;  winter  caae. 


ImmO  TARGET  ECHO  LEVEL 


TNO  report 


Pige 

33 
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Fig.  9;  Historical  sound  speed  profile  of 
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REACT,  using  ALMOST;  HMS;  dealer/  summer  case. 
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